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We have observed the orbital ordering in the ferromagnetic Mott-insulator LU2V2O7 by 
the polarized neutron diffraction technique. The orbital ordering pattern determined from 
the observed magnetic form factors can be explained in terms of a linear combination of 
wave functions \yz), \zx) and \xy); 

|0) = \[\\xy) + ^\yz) + \j\\zx) 

oc \{x + y + zf -r 2 ), 

where each orbital is extended toward the center-of-mass of the V tetrahedron. We discuss 
the stability of the ferromagnetic LU2V2O7, using a Hubbard Hamiltonian with these three 
orbitals. 
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1. Introduction 

Vanadium pyrochlores (RE)2V207, where RE = Lu, Yb and Tm, exhibit ferromagnetic 
and semiconducting states. 1 They crystallize in a face-centered-cubic structure with the space 
group Fd3m. The crystal is characterized as a three-dimensional network consisting of corner- 
sharing tetrahedra of V 4+ ions (Fig. 1). 

Lu 2 V 2 7 (a = 9.932 A) is a ferromagnetic Mott-insulator with Tc — 73 K. The electronic 
configuration of the V 4+ ions is fag) ■ The origin of the ferromagnetism still remains un- 
solved. 2 In order to understand the origin, we have to understand the electronic ground state 
of this system. 

A systematic distortion of VOg octahedra in a pyrochlore produces a trigonal crystal field 
to the vanadium ions 3 , suggesting that the ground state is the doubly degenerate e g orbital. 
This idea naturally leads to the existence of "anti-ferro" orbital ordering. In this article we 
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use the definition of orbital ordering in a broad sense. Namely, the orbital ordering means 
orbital polarization irrespective of its origin. On the other hand, band structure calculations 
showed that the non-degenerate a\ g has the lowest energy state, which is split from the t2 g 
orbitals in a cubic crystalline field. 2 

In this article, we present an orbital ordering pattern of LU2V2O7 obtained by the polarized 
neutron diffraction method. 

2. Experimental 

A single crystal of LU2V2O7 was grown by the floating-zone method in an Ar atmosphere. 
Starting materials, LU2O3 (4N) and V2O4 (3N), were mixed in a given molecular ratio and 
pressed into a rod. We obtained a single crystal with 4 mm in diameter and 10 cm in length. 
The samples used in the present experiment were cut to 4 mm in diameter and 2 mm in 
thickness. 

The polarized neutron diffraction measurements were performed using the Tohoku Univer- 
sity triple-axis spectrometer TOPAN in the JRR-3M reactor of JAERI at Tokai. Heusler alloy 
was used as a polarizer and energy of the incident neutrons was set at 80 meV (A = 1.011 A) 
or 70 meV (A = 1.081 A). The beam collimation was Open-100'-60'-Open. We also inserted a 
sapphire filter in front of the 2nd collimator to eliminate neutrons with higher incident energy. 

The data were taken at 3 K under an applied magnetic field of 3 T parallel to the (110) 
-axis of the single crystal specimen. Correction for the incomplete polarization (neutron po- 
larization P n ~ 94.7 % or polarization ratio R ~ 17.9) and estimations for the statistical un- 
certainty were properly carried out. 

We have performed similar measurements for two samples with different thickness and 
found that the extinction effect is negligible. Multiple scattering was also carefully checked by 
comparing the data taken at different incident neutron energy. 

3. Results 

In order to obtain the magnetic form factor, we used a well-known relationship in the 
polarized neutron diffraction method; the observed polarization ratio R which is the ratio of 
the diffracted intensities upon reversal of the incident neutron polarization direction, is related 
to 70 = Fm/Fn after the instrumental corrections as 



values were already obtained by the x-ray diffraction. 

Figure 1 shows the V 4+ ions in the unit cell of LU2V2O7. Since we have four sites in a 
tetrahedron, it is possible to assume that there are 4 types of V 4+ wave functions. Therefore, 
four types of magnetic form factors /1, /2, /a and f± are assumed. 




(1) 



where Fm and Fn are the corresponding magnetic and nuclear structure factors. The Fn 
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Using these magnetic form factors, magnetic structure factor can be expressed as 

j 

fi + h + h + U (all h, k, l = An) (a) 

fi - h - fs + h , , 

(h + k or k + l 

/1 + /2-/3-/4 > , L A . Q>) (2) 

or I + h = An) 

/1-/2 + /3-/4 J 

(all h,k,l = odd), (c) 

First, we consider the structure factors of the case (b). When there is no deviation from a 
spherical electric charge distribution (i.e., there is no orbital ordering), the magnetic structure 
factor Fm becomes because f\ = fi = fs = fi and no neutrons were scattered. Therefore, 
we can confirm the existence of the orbital ordering by observing the reflections at the Bragg 
points of the case (b). 

The results of the polarized neutron experiments are given in Table I. In this table, (///) m 
is the observed value, and (nf)o is the value after correcting the incomplete polarization. 

The final experimental results of the magnetic form factors are shown in Fig. 2 with 
open squares. The bars indicate the statistical errors. Note that the non zero values of the 
magnetic form factors of the (0 2 2), (4 2 2), (0 6 6), (8 2 2), (4 6 6) and (0 10 10) reflections, 
corresponding to the case (b) in eq.(2) and shown with arrows in Fig. 2, is a clear evidence 
for the existence of the orbital ordering in LU2V2O7. 

4. Calculation and Discussion 

We now proceed to calculate the magnetic form factors, which depend on the aspherical 
spin distribution, in order to compare the theoretically suggested orbital ordered configuration 
of the t2 g electrons with the experimental results. 

4-1 Magnetic Form Factor 

The electronic configuration of V 4+ ions is (t2g) 1 and the wave function of the ground 
state can be expressed in terms of a linear combination of the wave functions of t<i g - 

I*) = a\yz) + b\zx) + c\xy), (3) 

where \yz), \zx) and \xy) denote the wave functions of irreducible representation of t2 g . a, b 
and c correspond to the relative ratio of the wave functions, where a 2 + b 2 + c 2 = 1. 

The magnetic form factor f(K) is a quantity given by the Fourier transformation of the 
magnetic moment density distribution: 

f(K) = {*\* ■ Qjvp), (4) 

where er denotes an incident neutron spin, Qj_ = K x Q x K with K being the unit scattering 
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vector. The scattering operator Q is given by 

Q = e iK - r s + ^(lF + n), (5) 

where s and I are the spin and the orbital momentum operators for the electrons in V 4+ ions. 
T is a function defined as T = jij^pp J^ K r xe x dx as given in ref. 4. We ignore the second 
term in eq. (5) because the crystalline field may quench the orbital angular momentum and 
the magnetization of V 4+ ions. Therefore f(K) can be expressed as 

f(K) = (*\e* K - r s ± \*) 

= O'o)-^[(a 2 + 6 2 -2c 2 )(3cos 2 0-l) 

+ 12c sin 9 cos 9(a cos 4> + b sin 0) 

- 3sin6>{(a 2 - 6 2 )cos 2(f) - 2ab sin 20}] (j 2 ) 

- — [(4a 2 + 46 2 - c 2 ) (35 cos 4 6 - 30 cos 2 9 + 3) 
56 

+ 20c sin 9(7 cos 3 9 — 3 cos 9) (a cos <$> + 6sin0) 
+ 20sin 2 6»(7 cos 2 9-1) 

x {(a 2 - 6 2 ) cos 20- 2afrsin20} 
+ 140c sin 3 9 cos 9 (a cos 3<f> — b sin 3<f>) 

+ 35c 2 sin 4 6lcos40]O4), (6) 

where 9 and 4> denote the spherical coordinates of the scattering vector K relative to the 
quantization axes. (j n ) (n = 0, 2 and 4) are the Freeman- Watson radial functions, 5 which are 
calculated using the program cited in ref. 6. 

4-2 Orbital Ordered Structure in LU2V2O7 

To reproduce the magnetic form factor observed in LU2V2O7, we have made model calcu- 
lations on possible two types of models for orbital ordering structures. One is a model with 
three-fold degenerate orbitals in a tetragonal crystal field (model A) and the other is a model 
in a trigonal crystal field due to an oxygen lattice distortion (model B). 

As shown in Fig. 1, the unit cell of LU2V2O7 is constructed from a tetrahedron of four 
corner-shared V06 octahedra, and each V 4+ (3d 1 ) ion has a three- fold degenerate t 2g orbital in 
the tetragonal crystal field. The hopping amplitude of the electon depends on the direction of 
the hoppings and the orbital states due to the symmetries and anisotropics of the 3d orbitals. 
We define the local quantization axis on each V 4+ as shown in Fig. 3 (a) and express the t 2g 
orbital states on the site i as \xy)i, \yz)i, and \zx)i. Geometrically, the t 2g orbital on each site 
is approximately parallel to the surface of the V 4+ tetrahedron. Since the V 4+ ions structure 
is constructed with the corner-shared octahedra, the largest contribution for the t 2g electron 
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hopping is in general considered to be from t2 g -p-n-t2 g couplings. 7 Let us focus on the electron 
hoppings between the site 1 and 2 (Fig. 3 (a)). For example, the \yz)i and \zx}2 orbitals are 
roughly on the surface 123 in Fig. 3 (b). If the extents of the rotation and tilting were not very 
large, electron transfers between the orbitals parallel to the surface would be much larger than 
the other components. In fact, if the V-O-V angle were 9 = 180°, there would be a hopping 
process t\ from a t2 g -p-K~t2 g coupling between \yz)± and \zx)2- In the same way, the electron 
hopping from \zx)\ to \yz)2 is also favorable, because these orbitals are approximately parallel 
to the triangle 124 (see Fig. 3 (c)). On the other hand, since the direction from site 1 to 2 is 
nearly perpendicular to both the \xy)\ and \xy}2, the hopping amplitude from these orbitals 
should be small (see Fig. 3 (d)). In this way, there are two dominant hopping processes on each 
bond. Strictly speaking, the bond angle 9 in LU2V2O7 is smaller than 180° and there are two 
kinds of hopping processes between t2 g orbitals ti and t2 as shown in Fig. 4 (a). Nevertheless, 
since the angular dependence of the hoppings are ti oc cos 9 — 1 and £2 oc cos 9 + 1, we can 
neglect the hopping process £2 as a first order approximation. 

Let us consider the model A (see Fig. 4 (a)) based on the Hubbard Hamiltonian with 
three-fold degenerate orbitals in the limit t±/U <C 1, where we only take into account the 
most dominant hopping process t±. Since the electron in the \yz)i orbital can move only 
within the hatched surface in Fig. 3 (a), i.e., the triangle 123 of the tetrahedron, only one 
electron can exist in the three orbital states \yz)±, \zx)2, and \xy)$ on the surface-plane to gain 
the kinetic exchange energy. Thus due to the Pauli principle it is expected that the energy gain 
is maximum when the 3d electrons on V 4+ ions occupy the t2 g orbitals which lie on different 
surface-planes so that the virtual hoppings of electrons t\ are not prohibited. An example 
of such an orbital ordering is realized when all occupied orbital states on the V 4+ ions can 
simultaneously be expressed as one orbital, e.g., \yz) orbital. Namely a "ferro" -orbital ordered 
state should be favored. Note that the orbitals on each site face the different directions even 
in the "ferro" -orbital ordered state since we define the local quantization axis on each site 
differently. 

Next, we consider the model B (see Fig. 4 (b)). In the trigonal crystal field, the three-fold 
degeneracy of t2 g orbitals is lifted. The lowest energy state which has an a\ g symmetry is 
written as, 

|0> = ±={\xv) + \yz) + \zx)). (7) 

The orbital ordering structure is shown in Fig. 5, where each orbital is extended toward the 
center-of mass of the tetrahedron. 

Now we compare the model A and B. It is expected that there is a competition between 
kinetic energy gain through the orbital exchange and splitting energy gain due to the trigonal 
crystal field, which determine the orbital structure in LU2V2O7. Assuming the orbital states in 
the models, we calculate the magnetic form factors and compare the results with that observed 
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in the experiment. As shown in Fig. 2, the form factors in the model B agree much better with 
the experimental results. The method of least squares indicates that "ferro" -orbital ordered 
state with 

\4>)i = a\xy)i + b\yz)i + c\zx)i, (8) 

where a = 0.67 ±0.25, b = 0.40 ±0.23, c = 0.62 ±0.14, is favored. The orbital state \4>}i agrees 
with that observed in model B (a = b = c = 0.58) within the computed error bars. Thus 
we conclude that the orbital ordered structure shown in Fig. 5 is realized in LU2V2O7. Our 
orbital structure model is also consistent with that proposed by Shamoto et al. In this way, 
it is expected that the effects of the trigonal field is bigger than that of the exchange energy 
of the electrons in LU2V2O7. In fact, the temperature T OQ at which orbital ordering occurs is 
much higher than the Curie temperature T c = 73K. 

Finally we check the stability of the ferromagnetic state, which has been observed in the 
experiment, on the model B. Following the procedure in Ref. 8, we start from the Hubbard 
Hamiltonian with 3d orbitals: 

f/ - V V f m ' m V f r , 
i,j m,m',(T 

± U ^ ] IT'iml T^iml ^ ^ ^ ^ima^im' — a 

+ (U'-J H ) 

— Jh ^2 ( C im y \ Cim l C im l j c *m'T 

+ AmA m i C im'\Cim'i ± H.C.) 

+ ^ K n n im(J . (9) 

i,m,cr 

Here U is an intra-orbital Coulomb interaction, U' an inter-orbital Coulomb interaction, Jh 
Hund's rule coupling, and A m splitting energy due to the trigonal crystal field. Indices m, m! 
indicate orbital states, i,j sites, and a spin states. 

Under the trigonal crystal field, the three-fold degeneracy of the t2 g orbital is lifted: the 
lowest energy state |0) which is defined in eq.(7) and the excited states with the e g symmetry. 
We define that the splitting energy for the lowest energy state |0) is zero, and that for the 
excited state Ao in the Hamiltonian (9). Given the transfers t[ and t' 2 as shown in Fig 4 (b), 
we can obtain an effective Hamiltonian in the second order perturbation for t'i/U and t' 2 /U, 

I6tf 8t' 2 2 ( J H \ \ fl 



Heff = ~ 1 ~~9U + w 1 V 1 T7^~ J J V 4 s i " s 2 
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) 



(10) 



where we assume the condition that U, U' 3> Jh 3> A = 0. When the spins of the electrons 
on site 1 and site 2 are antiparallel, the energy gain through the virtual hopping process to 
the intermediate states |0) and |±) are reduced by the first term of eq.(10). On the other 
hand, when the spins are parallel, the transfer t\ is forbidden by Pauli's principle and the 
energy gain through the intermediate states |±) are possible by the transfers t' 2 as given in 
the second term of eq.(10). We have checked the stability of the ferromagnetic state, assuming 
that U = 6.0 eV. 2 The result indicates that the ferromagnetic state is stable for Jh ^ 1-8 
eV, where we use the relation \t' 2 \ = ji'J/2 (see Appendix). Hund's rule couplings estimated 
in Lu 2 V 2 7 J H = 1.5 eV (U' = 3.0) 2 and J H = 0.7 eV (W = 4.6 eV) 9 are slightly smaller 
than the critical value 1.8 eV. This may contradict the experimental results. 

So far we take into account the hopping processes between ti g orbitals. However, including 
the effects of the hopping process from t2 g orbitals to e g orbitals (\x 2 —y 2 ) and \3z 2 —r 2 )) might 
be important to stabilize the ferromagnetism in t 2 g orbital system (see fig. 4(b)), as predicted 
by Mochizuki and Imada in the Perovskite-type Ti oxides. 10 In fact, even in LU2V2O7, the 
ratio of the hopping processes for the V-O-V angle 9 = 130° are |i 3 |/|i]J ~ 3 (see Appendix). 
Therefore the effects of the virtual hopping process t' 3 can not be neglected even if there is a 
level splitting Ai due to the crystal field. Including the effects of i 3 , the second order effective 
Hamiltonian is written as 



J( Jh, 0, Ao, Ai) has been calculated as a function of Jh and Ai with fixed values 8 = 130° 
and Ao = 0.4 eV. The parameter ranges where the ferromagnetic state is stable are shown in 
Fig. 6. In this way, the energy gain through the virtual hopping process t 3 can stabilize the 
ferromagnetism rather than the antiferromagnetism. It is natural to think that in LU2V2O7 
the virtual hopping process i 3 can play an important role to stabilize the ferromagnetism, 
since the ferromagnetic ground state has been observed experimentally in this material. In 
fact, assuming the level splitting energy Ai = 2.0 eV, 9 both estimated values of Hund's rule 
coupling Jh = 1.5 eV 2 and Jh = 0.7 eV 9 indicate the presence of the ferromagnetic ground 
state in LU2V2O7, which is consistent with the experimental results. Once the V-O-V bond 
angle is away from 180°, the transfer amplitude between ti g to e g orbitals can not be negligible. 
As a result, the ferromagnetic state can be favorable due to the energy gain by the second 
order perturbation processes through the ti and £3 hopping transfers. 
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Appendix 

Under the tetragonal crystal field, a five- fold degeneracy of 3d orbitals on V 4+ (3d 1 ) ion 
is lifted to three- fold lower levels t2 g state and two- fold higher levels e g state (see Fig. 4 (a)). 
Because of symmetries and anisotropics of the 3d orbitals, the hopping amplitude depends 
on the direction of the hoppings and the orbital states. Let us consider the hopping process 
from \yz)i to the orbitals on the site 2. Since the V-O-V bond angle in LU2V2O7 is about 
9 = 130°, the \yz)i orbital hybridizes with not only \zx}2 but also \yz}2 and \3z 2 — r 2 )i- Thus 
hopping amplitudes £2 and £3 defined in Fig. 4 (a) are non-zero. The angle dependences of the 
hoppings are written as 

h = -^V pd7T V pd7T {cos6- 1), (12) 
t2 = -^V pd7T V pd7T (cos 9 + 1), (13) 

h = -—^Vpd^Vpda sm6, (14) 

where A is the charge-transfer energy which describes the energy difference between occupied 
O 2p and unoccupied V 3d levels, and V pda and V pdn are interatomic matrix elements for a 
and 7r bonds, respectively. 

Under the trigonal crystal field, the degeneracy of ti g states is lifted. The lowest state |0) 
is defined in eq.(7) and two states with the e g symmetry |+) and |— ) are 

|+) = "^(N> + u\yz) + uj 2 \zx}), (15) 

\~) = -^\xy) +^ 2 \yz) +uj\zx)), (16) 

where to = exp(27ri/3). Note that the higher e g orbitals \x 2 — y 2 ) and \3z 2 — r 2 ) remain unlifted. 
Using the hopping processes ti, £2, and £3, possible hopping processes shown in Fig. 4 (b) have 
been calculated and the results are 

2 

t[ = —^V pd7r V pd7r cos 9, (17) 

t>2 = 3~A Vpdn Vpdn C ° S0, ^ 

4 = -^^VpdvVpdv sm9. (19) 

The ratio of interatomic matrix element V pda and V pdn is fixed at V pda /V pdn = —2.15, which 
is the value obtained in the Perovskite compounds. 11 Using the V-O-V angle 9 = 130°, the 
ratio of the hopping processes are |^| : \t' 2 \ : \t' 3 \ ~ 2 : 1 : 6. 
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Figure.l V 4+ ions in the unit cell of (RE) 2 V 2 7 (RE = Lu, Yb and Tm). 

Figure. 2 Comparison between the observed and the calculated magnetic form factors. The 
open squares with bars correspond to the experimental values. The open circles are the 
calculated values of the orbital ordering model A, and the filled circles are the ones of the 
model B. The computed error bars of the least squares calculation are within plots. 

Figure. 3 (a) The definition of the directions of the x, y, z axes at each V site on a tetrahedron. 
Circles represent O ions, (b) Orbital states on the site 1 and 2, which are approximately 
parallel to the hatched plane 123. The i 2s electron hoppings through t2 g -p n -t2g couplings 
are shown by arrows (p orbitals are omitted for simplicity), (c) Orbital states on the site 
1 and 2, which are approximately parallel to the plane 124. The t 2g electron hoppings 
are shown by arrows, (d) An orbital state on the site 1 (2) which is approximately par- 
allel to the plane 134 (234). The i 2a electron in these orbitals can not hop from site 1 to site 2. 

Figure. 4 Characteristic hopping process, (a) in a tetragonal crystal field, (b) in a trigonal 
crystal field. 

Figure. 5 Schematic view of the orbital ordering in Lu 2 V 2 07 (model B). 

Figure. 6 The parameter region where ferromagnetic (antiferromagnetic) interaction couplings 
are realized. Estimated values by Shamoto et al. and Mizokawa are shown by circles. 
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(a) in a tetragonal crystal field (model A) 




(b) in a trigonal crystal field (model B) 



Fig. 4. 
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Fig. 5. 



H. Ichikawa et. al. 



15/17 



J. Phys. Soc. Jpn. 



Full Paper 



- 


FM 


i 

Shamoto et al. 
• 


- 




Mizokawa 






AF 

1 



1.5 2 2.5 

4x(eV) 

Fig. 6. 



H. Ichikawa et. al. 



16/17 



J. Phys. Soc. Jpn. 



Full Paper 



Table I. Comparison between the experimental and calculated data. 

(pf)m are the observed values, (/i/)o are the observed values after correction, (/i/)^, (fJ,f)s are the 
calculated values, corresponding to the model A and B with /i = 0.63 fiB- Parenthesis are the statistical 
errors. 
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